Rayleigh superradiance and dynamic Bragg gratings in an end-pumped Bose-Einstein condensate 
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We study experimentally superradiant Rayleigh scattering from a Bose-Einstein condensate (BEC) in a new 
parameter regime where pump depletion and the exchange of photons between the endfire modes are important. 
Through experiments and simulations we show that collective atom light coupling leads to the self-organized 
formation of dynamic Bragg gratings within the sample. These gratings lead to an efficient back-scattering of 
pump photons and optical resonator structures within the BEC. 
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With its extremely high optical depth and unique coher- 
ence properties, a Bose Einstein condensate provides an ideal 
object on which to study collective light scattering with the 
goal of generating and probing light-atom correlations. Su- 
perradiant light scattering (SLS) from ultra-cold atomic en- 
sembles has been recognized as a method to generate entan- 
gled atoms and photons due to the fact that the interaction 
Hamiltonian has the generic form of a parametric amplifier 
H cx ivb', where correlations based on momentum conserva- 
tion arise in the case of Rayleigh scattering [1], and angular 
momentum in the case of Raman processes yfl. Entanglement 
in the ultra-low gain regime of such an interaction, when far 
less than one atom-photon excitation pair is generated on av- 
erage, forms the basis for a quantum repeater [3] and has been 
studied extensively Q4|] . Here, we are interested in the high 
gain - superradiant - regime where the detection of entangle- 
ment requires in general the measurement of both the phase 
and amplitude of the light and matter waves at the sub-shot 
noise level. While this is typically hampered by the lack of 
an atomic analogue to homodyne detection of light, a recent 
proposal claims entanglement may be detected by just count- 
ing photons and atoms |5[]. Nonetheless, the full dynamics 
of coupled matter and light waves where one does not place 
limitations on the depletion of either the condensate or input 
light, and where propagation effects are considered, remain 
important issues if such correlations and entanglement are to 
be made useful resources. 

To this end, we explore superradiant Rayleigh scattering 
from a trapped, cigar shaped BEC as the pump detuning is 
varied while the single particle scattering rate R is kept con- 
stant. In this way, we investigate the effect of the detun- 
ing of the pump beam in the process, and move between the 
case where the pump beam remains essentially undepleted by 
the scattering, to the situation where superradiant scattering 
is 'clamped' by a lack of photons in the pump beam. Cru- 
cial to these dynamics is the structure that builds up along 
the long axis of the condensate, demonstrating characteristics 
from 'Dicke' superradiance from extended samples a 
fact reco gniz ed recently both experimentally II9L llOfl and theo- 
retically HHd^dll 0]. Contrary to most earlier experimen- 
tal work on the subject J<| [l5l Q3, \Ym . we end-pump along 
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the long axis with a beam mode-matched to the transverse 
cross-section of the BEC in order to optimize the coupling 
between pump photons and atoms. As opposed to the study in 
01811 . which also probed the BEC along its long axis, we use 
the back-scattered light to monitor the superradiant dynamics 
with very good time resolution. In view of our ultimate goal 
to probe light-atom correlations, we have performed a care- 
ful calibration of the detection system. As a consequence, we 
find quantitative agreement with simulations of a ID semi- 
classical model over a wide parameter range. The quality of 
the agreement with experiment supports the use of the simula- 
tions to understand the spatial and temporal dynamics inside 
the sample. 

Like superradiance (SR) in electronically inverted samples 
lfl9ll . superradiant light scattering from an ultra-cold atomic 
sample is a process whereby an initially unoccupied elec- 
tromagnetic field mode becomes weakly populated through 
spontaneous emission; these modes are then amplified with 
the highest gain along the direction of greatest optical depth. 
In the end-pumped geometry (see Fig. [TJ, photons are ab- 
sorbed from the forward travelling beam E + and scattered 
into the backward travelling beam £L, leading to a concomi- 
tant change in the scattered atom's momentum of 2Tiki and 
a recoil shift in the back-scattered light of four times the re- 
coil frequency u) r . The forward and backward travelling light 
waves interfere to give a moving intensity grating with a spa- 
tial period of A/2 that varies in visibility and phase over the 
length of the sample. In the case of Rayleigh scattering, the 
internal ground state of the atom does not change in the pro- 
cess, leading to a density modulation due to interference be- 
tween the different momentum orders that has the same period 
as the light intensity modulation. Thus, the physical picture is 
one of four wave-mixing. 

In our experiments, superradiant Rayleigh scattering was 
induced in a trapped BEC by illuminating it with a pulse of 
off-resonant light along the long axis of the condensate, as 
shown in Fig. Q] The BEC was generated by evaporatively 
cooling a cloud of ^^Rb atoms in the \F = l,m F = — 1) 
hyperfine state in a Ioffe-Pritchard magnetic trap. We obtain 
cigar shaped condensates containing 1.35 x 10 6 atoms, with 
in-trap Thomas Fermi radii of r± = 6.4 and r|| = 65 /xm in 
the radial and axial directions and with no discernible thermal 
fraction. The pump light was detuned by a variable amount 
from the \F = l,m F = -1) -> \F = 2,m F = -2) 
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FIG. 1: (Color online) (a) The BEC was end-pumped by a focused, 
circularly polarized laser beam, and the back-scattered light was 
directed through a quarter wave plate (L/4) and polarizing beam- 
splitter (PBS) onto a sensitive photodetector (PD). (b) Traces from 
PD are shown for low and high pump powers at a detuning of 
S = — 2n x 2.6 GHz. Simulations for the same parameters are pre- 
sented as red dash-dot lines. Absorption images after 45 ms time of 
flight (TOF) of the corresponding atomic distributions are shown in 

(c) (low power, pump pulse duration 200 /is) and (d) (high power, 
pump pulse duration 50 /is). Circles indicating the separation of ad- 
jacent momentum orders after 45 ms TOF are shown in (d); note that 
they originate from the input light end (x = 0) of the BEC. 

transition on the Dl line of ^Rb at 795 nm, and circularly 
polarized with respect to the long axis of the trap. All data 
presented is for red detunings (6 — oji — ojq < 0); rectan- 
gular pump pulse envelopes; and where the atoms were inter- 
rogated in-trap, with the trapping potential extinguished im- 
mediately after the end of the pump pulse. The beam was 
focused to a waist of 13 /im at the center of the condensate 
with negligible change of beam size over the length of the 
BEC. Light was back-scattered by the sample in the same po- 
larization as the input beam, and thus the backward travel- 
ling light was reflected by the polarizing beamsplitter, then 
directed onto a sensitive PIN diode photodetector. The detec- 
tor, with a bandwidth of 400 kHz, is shot-noise limited for 
photon fluxes greater than 10 5 /is -1 . To avoid back reflec- 
tions from optics and cell windows that seed the process, the 
pump beam was inserted at a slight angle (less than 2°). Pic- 
tures of the atoms were obtained after 45 ms time of flight by 
resonant absorption imaging. FigureQTb) shows experimental 
and simulated time traces for high and low pump powers, with 
the corresponding atomic distributions shown in Fig. [Tfc) and 

(d) . For single particle scattering rates much smaller than the 
recoil frequency uj r = 27r x 3.6 kHz, scattering to higher 
atomic momentum orders occurs sequentially on a time scale 
~ r r = 2tt /uo r . Figure [Tfc) shows this case, where the trans- 
fer is limited to the first order. When R ~ u> r , there is suf- 
ficient gain for atoms to be back-scattered into negative mo- 
mentum orders, i.e., the Kaptiza-Dirac regime where atoms 
absorb back-scattered light and re-emit into the forward di- 
rection, as is evident in Fig. Old). An asymmetry in distance 
between forward and backward scattered atoms and the cen- 
ter of the original condensate after time of flight, observed 
previously in Q9J, [TT1], is visible in Fig.Q~fd). This asymmetry 
can be traced back to the spatial inhomogeneity of superradi- 



ant scattering favoring the input end of the condensate, where 
the amplitude of the reflected light E- is highest. The spatial 
dimensions of the condensate and the slow expansion upon re- 
lease from the trap along the long axis are such that this spatial 
feature of the scattering is evident after 45 ms time of flight. 
The momentum distributions are somewhat distorted due to 
the input angle of the beam: when the beam is aligned paral- 
lel to the long axis of the BEC, the Bragg condition is satisfied 
for the atoms centered around zero transverse momentum, but 
at the slight incident angle used, the patterns become more 
complicated. Nevertheless, the arrival times and amplitudes 
of the first superradiant light pulses are very well described 
by the simulations. 

The dynamics are simulated by numerically solving ID 
Maxwell-Schrodinger equations for the evolution of the 
atomic wavefunctions and the incident and generated light 
fields [6, 20]. The situation of SLS from a BEC was derived 
in ifTTl 12], and we follow closely their notation. In the end- 
pumped geometry, the Maxwell-Schrodinger equations read: 
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such that ip(x,t) = £m V'mOM)e~ 4( ''' mt ~ m ' C!:E) is the 
slowly varying ID wavefunction of momentum or- 
der m — 2n for integer n, with cj m = m 2 oj r ; and 
the electric field (suppressing polarization) is given by 
E = E + (x,t)e- l(uit ' k ' x ^ + E-{x,t)e-^ lt+klX *> + /i.e.; 
where E±(x,t) are slowly varying envelopes. The light 
fields are scaled such that E± = e±yjTiui r ki/ £qA\ and the 
scaled space and time variables are related to real space 
parameters by £ = kix and t = 2u r t. The coupling 
constant describing the strength of the interaction is given 
by A = 1/4 • T/S ■ <jq /A, where V denotes the linewidth of 
the optical transition, cto the resonant atomic absorption cross 
section, and A the cross section of the interaction region. 
Retardation effects have been neglected in Eqns. (|2]) and (0) 
given the length of the condensate 2ry = 130 /im, which 
allows us to discard a time derivative term. In this approxi- 
mation, these equations can be considered as describing the 
self-consistent light field distribution in a continuous array 
of Bragg gratings formed by the density modulation of the 
matter wave. In this ID model, the endfire mode described 
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by e+ and the pump mode coincide. A ID treatment of the 
problem is supported by the fact that the Fresnel number 
of the condensate is of order one, which suggests that the 
emission of light is predominantly into a single transverse 
mode [211]. The harmonic trapping potential and the mean 
field interaction term describing collisions between atoms 
have been omitted, because they have a limited effect on the 
focus of these experiments - the first superradiant pulse [22]. 

Equations (0, and describe a Raman interaction 
where a ladder of momentum states is coupled by two counter- 
propagating light fields. It is terms three and four that domi- 
nate in Eqn. ([TJ): these describe the local coupling to the near- 
est momentum states via exchange of photons between e + and 
e_. The final terms in Eqn. <[TJ account for phase rotation 
of the matter wave due to the light shift. Equations 01 and 
01 show terms equivalent to the coupling terms in Eqn. (Q}. 
Specifically, the growth of e_ occurs with the corresponding 
growth of recoiling atoms VVn+2 and the decrease of e + pho- 
tons and tp m atoms. The last terms in Eqns. and de- 
scribe the effect on the light of the slowly varying refractive 
index due to the large scale atomic density distribution. 

To compare experimental results with simulations, the 
equations (0, (0 and (0 were solved numerically for exper- 
imental parameters. The initial wavefunction i[jq was taken 
to be a ID Thomas-Fermi profile normalized to the num- 
ber of atoms in the trap iV at . The boundary conditions for 
the light fields were typically taken to be e + (0, r) = ei 
and e_(fc;L,r) = 0. Note that |e_(0, r)| is proportional 
to the back-scattered light intensity measured in experiment, 
and that the total photon flux is conserved: |e+(0, r)| 2 = 
|e_(0,r)| 2 + \e + (kiL,r)\ 2 . As Eqns. (0, (0 and contain 
no explicit noise term to instigate superradiant scattering, we 
seed the process by taking a non-zero first order momentum 
component -02 = 4>o / V~Nau corresponding to a single delo- 
calized atom in the first side-mode [12]. The random nature 
of the initiation of superradiant scattering may be modeled in 
the present formalism by using random initial conditions but 
here we focus on the large scale dynamics of the system rather 
than its noise properties 12011 . 

Figure 0a) shows, for both experimental data (symbols) 
and simulations (lines), the peak values of the first superradi- 
ant pulse as the detuning of the pump beam was varied while 
keeping the single particle scattering rate R constant. Re- 
sults are shown for four values of R in the range R <C uJ r 
to R ~ u> r . For a large portion of the graph, the peak value is 
essentially independent of the detuning, and there is excellent 
agreement between simulations and data. The input powers 
in the simulations have been scaled up by a common factor of 
22% B23I1 . For lower detunings, the peak value of the emitted 
pulse falls away, and the experimental data reaches our detec- 
tion resolution for 6 « — 2tt x 300 MHz, while the simulations 
show the same qualitative behaviour. 

To further study the two regimes, Fig. 0b) shows experi- 
mental time traces for several different values of R at a high 
detuning, and Fig. 0c) shows experimental traces for four 
low detunings while R is kept constant. Figure 0b) confirms 
the qualitative features evident in Fig. 0b): in the high detun- 
ing regime, the superradiant peaks arrive earlier and are more 
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FIG. 2: (Color online) (a) Experimental (symbols) and simulated 
(lines) peak photon flux of the first superradiant pulse as a func- 
tion of detuning for R —2.2 (circles, solid line), 6.4 (right trian- 
gles, dashed line), 10.7 (upward triangles, dotted line), 15.9 (squares, 
dash-dot line)xl0 3 s _1 . Errorbar limits are the standard error of 
the mean of 5 realizations, and for clarity are omitted in the low 
detuning portion of the graph, (b) Traces corresponding to the flat 
portion of Fig. |2fa). The detuning was S = — 2tv x 2.6 GHz, for 
R = 2.2, 4.3, 6.4, 8.6, 10.7, 15.9 x 10 3 s~\ For clarity, the traces 
are edited to show only the first pulse, (c) Traces corresponding to 
the decaying portion in Fig. [2a). For 7? = 2.2 x 10 3 s _1 and detun- 
ings S = -2ttx 2.0 (black), 1.0 (red), 0.7 (blue), 0.5 (green) GHz. 

sharply peaked the higher the pump power. Simple models of 
superradiance predict a pulse height proportional to R and a 
pulse delay inversely proportional to R O], while we observe 
a faster increase of the peak power and a slower decrease in 
pulse delay both experimentally and in the simulations. Fig- 
ure0c) shows how the superradiant pulse intensity is clamped 
at low detunings. The lower amplitudes can be traced back to 
the lower incident photon fluxes and increasingly important 
pump depletion in the low detuning case. The transition from 
high detuning behaviour to the pump-depletion regime occurs 
in the experimental data at points where the incident photon 
flux is approximately 10 times the peak reflected photon flux. 
Above this point, while the amplitude of the reflected pulses 
drops, the observed peak reflectivity of the sample increases 
sharply, up to values of ~ 30%. The simulations show the 
same qualitative behaviour in this regime, and we attribute the 
loss of quantitative agreement to the incoherent losses that are 
not captured by the model: i.e., emission into different modes. 
The spatial dependence of the light and matter waves plays a 
critical role here, and to gain more insight into the behavior 
of the system, we explore the light and matter wave dynamics 
inside the sample through simulations. 

Figure [3] shows the results of simulations for the lowest 
pump power case for high (left side) and low (right side) 



4 




24 L= , ^ , 1 I ■ — 

50 100 50 100 

Position (|im) Position (um) 



FIG. 3: (Color online) Snapshots of system dynamics inside the 
sample for high detuning 8 = — 2nx 2 GHz ((a) and (c)) and 
low detuning S = — 2tyx 294 MHz ((b) and (d)) at the same rate 
R — 2.2 x 10 3 s _1 after an interaction time t = 87 /is. (a) and 
(b) show atomic densities in the momentum modes: \1p-2 ^(dash- 
dot line), |i/)o| 2 (solidline), |02 | 2 (dashed line), 04 \ 2 (dotted line), (c) 
and (d) display light flux (oc \e± | 2 ) within the sample for pump light 
e+ (solid lines; in (c), red line and left axis) and reflected light e_ 
(dashed lines; right scale in (c)). Note the shared scale of e_ in (c) 
and both e+ and e_ in (d). 

pump detuning, after 87 of interaction time with the pump 
beam. Along the long axis of the BEC, atomic distributions 
are shown in the upper row ((a) and (b)), and the scaled in- 
tensities oc |e+| and |e_| in the lower row ((c) and (d)). 
The general dynamics for low input power, high detuning are 
as follows: The back-scattered light intensity in the sample 
builds up at the input end because there it sees gain from ap- 
proximately the entire length of the BEC. At this point, atoms 
are transferred from 0o to 02 concurrently with the growth of 
e_ and reduction in e+. When the population in 0rj is suf- 
ficiently depleted at the input edge of the BEC, the process 
slows down, and the light field envelopes move towards the 
centre of the condensate, where |0o| is still large, and the ex- 
change between the two light fields continues; this is the time 
shown in Fig. [3] At this time, |^o| grows again at the input 
end of the condensate, driven there by the destruction of e_ 
photons generated further inside the sample, and 02 atoms. In 
this way, the back-scattered photon flux out the input end of 
the condensate stops, and the first superradiant pulse has been 
emitted. The basis of 'ringing' behaviour - the emission of 
subsequent SR pulses - is merely a repetition of the dynam- 



ics described above. Furthermore, a fascinating implication 
of the above dynamics, visible in Fig.[3](c) and (d), is that at 
times the light intensity within the BEC is higher than outside 
- the interaction leads to the formation of an optical resonator, 
where partially reflecting mirrors are formed by the density 
modulation due to the interference of stationary and recoiling 
matter-waves. These Bragg gratings are centered where 0o 
and -02 cross. 

The dynamics in the low detuning case is similar, but with 
two significant differences. Due to the increased light scatter- 
ing cross-section at low detunings, the pump light is signifi- 
cantly depleted in its passage through the BEC, as is evident 
in Fig. [3d). Accordingly, the build up of back-scattered light 
is even more localized at the input end of the BEC, and the 
scattering of 0o to 02 atoms is limited to a very narrow region 
of atoms. This narrowness is reflected at later times in the 
length of the 'resonator', as is evident in Fig. 0b). Thus, the 
spatial dependence of the pump depletion and hence the size 
of the scattering region leads to the reduced amplitude of the 
back-scattered pulses in Fig. [2] The second main difference is 
that the 04 mode becomes significantly populated at low de- 
tunings, as can be seen in Fig.[3|b). This means that the first 
SR pulse is not stopped by a lack of 0o atoms, but rather can 
continue due to the scattering from 02 to 04 via the destruc- 
tion of a pump photon. This can be seen in Fig. [3jd), where 
there is more light 'leaking' from the resonator compared to 
Fig. 0c); and in Fig. [2jc), where the low detuning pulses are 
broader, and show a secondary peak soon after the first. In 
general, the dynamics become complicated as the number of 
significantly atomic orders and hence number of timescales in 
the problem (cx l/w m ) increases. 

In summary, we have studied superradiant scattering of 
light by Bose condensed atomic samples in new parameter 
regimes and find very good agreement between experimen- 
tally detected light pulses and simulations. The simulations 
take into account the build-up of longitudinal spatial structure 
inside the sample of both the light and matter waves. The 
results help to identify a suitable parameter regime for the ex- 
perimental study of correlations between light and atoms that 
arise from superradiance. 
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